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Abstract 

We present the evidence of the low defect density at Ge/Ge02 interfaces in terms of first- 
principles total energy calculations. The energy advantages of the atom emission from the Ge/GeO"2 
interface to release the stress due to the lattice mismatch are compared with those from the Si/Si02 
interface. The energy advantages of the Ge/Ge02 are found to be smaller than those of the Si/Si02 
because of the high flexibility of the bonding networks in Ge02- Thus, the suppression of the 
Ge-atom emission during the oxidation process leads to the improved electrical properties of the 
Ge/Ge02 interfaces. 

PACS numbers: 81.65.Mq, 68.35.-p, 71.15.Mb, 81.05.Cy 
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For the past several decades, the performance of Si/SiC>2 based transistor has been mainly 
accomplished by scaling the transistor dimensions resulting in a higher circuit speed, higher 
packing densities, and less power consumption. At present, the aggressive scaling of the 
SiC>2 gate dielectric in silicon metal-oxide-semiconductor (MOS) devices leads to unaccept- 
ably high gate leakage currents. Thus, the scaling of advanced devices is approaching its 
technological and fundamental limits. To further improve the performance of complemen- 
tary MOS devices, alternative channel materials are considered. Germanium is introduced 
as one of the promising candidates for beyond scaling devices, because the intrinsic carrier 
mobility of germanium is higher than that of silicon. To realize high-performance devices 
with germanium channel, one of the most crucial issues is the formation of gate stacks with 
superior interface properties. Ge/Ge02 interface is one of the most important issue because 
it exists even in Ge/high-A; oxide interfaces. Although one might expect that germanium, as 
another elemental column- IV semiconductor, would behave similarly to silicon, the Ge/Ge02 
interface had been generally considered to be more defective than the Si/Si02 interface, in 
which the interface trap density typically lie in the range of the latter half of 10 11 ~ 10 12 
cm~ 2 eV _1 .1]. Very recently, Matsubara et al. 2| have reported that the minimum value 
of the interface trap density lower than 10 11 cm _2 eV _1 can be obtained for Ge/Ge02 MOS 
interfaces fabricated by dry oxidation without any hydrogen passivation treatment. In addi- 
tion, the current authors' group (T. H. and H. W.) also accomplished high-quality Ge/Ge02 
interfaces by conventional thermal treatment of germanium substrates 3|. 

On the theoretical side, the structural and electronic properties of the Ge/Ge02 interface 
have been investigated 0, O, Ej]. Houssa et al. [(J simulated the density of germanium 
dangling bonds at the Ge/Ge02 interface as a function of the oxidation temperature by 
combining viscoelastic data of Ge02 and the modified Maxwell's model, and claimed that 
the density of germanium dangling bonds is less than that of silicon dangling bonds. Their 
results are in the good agreement with the experiments done by Matsubara et al. [2] and our 
result The future development of the passivation technologies to deactivate the interface 
defects, by introducing terminators such as hydrogen, sulfur, silicon, and so on, will make it 
possible to achieve the lower interface defect density in Ge/Ge02- However, the formation of 
the defects in the Ge/Ge02 interface region during the oxidation process is not well-known. 

In this Letter, we implement the first-principles total energy calculation of the Ge/Ge02 
interface during the oxidation process. In the case of the Si/Si02 interface, Kageshima 
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FIG. 1: Atomic configurations of surface models, (a) No germanium atom is emitted from the 
structure which includes three oxygen atoms, (b) One germanium atom per unit cell is emitted 
from the structure which includes three oxygen atoms, (c) One germanium atom per unit cell is 
emitted from the structure which includes six oxygen atoms. The empty circles, the grey circles, 
and the filled circles are germanium atoms, hydrogen atoms, and oxygen atoms, respectively, and 
the broken circles indicate the position where the germanium atom is emitted. 

and Shiraishi |7|] simulated oxidation processes by first principles. They found that a quartz 
structure can be obtained if silicon atoms are kicked out from the interface during oxidation, 
and these emitted silicon atoms might result in the creation of the interface defects. We 
compared the emission probabilities of the atom from the Ge/GeC>2 and Si/Si02 interfaces in 
terms of the total energy calculation following the oxidation process proposed by Kageshima 
and Shiraishi. To the best of our knowledge, this is the first theoretical investigation of the 
evidence of the low interface defect density at the Ge/Ge02 interfaces by first-principles 
total energy calculations. We found that the energy advantages of the germanium-atom 
emission are lower than that of the silicon-atom emission, and the small energy advantages 
are strongly related to the high flexibility of the bonding networks in Ge02- 
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real-space finite-difference method [id, M, Il2l |. In this method, basis sets used to expand 



Our calculations were performed based on the density functional theory 

aa 

using; the 



wave functions are not required and boundary conditions are not constrained to be periodic. 
Our models are imposed the periodic boundary condition in the directions parallel to the 
surface and the isolated boundary condition in the direction per pen dicular to the surface. 
The norm- conserving pseudopotentials of Troullier and Martins \l<L Jjj] in the Kleinman- 
Bylander representation [l5] were employed to describe the electron-ion interaction and the 
exchange- correlation potential was treated within the local density approximation . The 
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FIG. 2: Atomic configurations of interface models in which additional three oxygen atoms are 
inserted into the Ge/Ge02 interface region (a) before and (b) after emission. The meanings of the 
symbols are the same with those in Fig. [TJ 

cutoff energy was set at 112 Ry, which corresponds to a grid spacing of 0.30 a.u., and 
a further higher cutoff energy was set at 1011 Ry in the vicinity of the nuclei with the 



augmentation of the double-grid technique 
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12]. Eight k-points in the lxl lateral 



unit cell were used for the Brillouin zone sampling. Figure [T] shows an example of the 
Ge(100) surface models. The initial structure of the surface model has nine germanium 
atomic layers, and the size of the supercell was L x = L y = and L z = 5.5ao, where 
L x , L y , and L z are the lengths of the supercell in the x, y, and z directions, respectively, 



171 ] . Both sides of the surface are 



and cio is the experimental lattice constant of the bulk 
simply terminated by hydrogen atoms. We sequentially inserted oxygen atoms between 
Ge-Ge bonds from the surface, assuming atomical layer-by-layer oxide growth, and finally 
introduced six oxygen atoms per unit cell. As the further oxidation model, Fig. [2] shows an 
example of the interface models. The initial structure of the interface model has three Ge02 
molecule layers and seven germanium atomic layers. The length of the supercell in the z 
direction was L z = 7.1ao, and the other computational details were the same with those 
described in the surface models. We calculated the interface models inserted three or six 
oxygen atoms per unit cell into the interface region. During the first-principles structural 
optimization, we relaxed all the atoms except the germanium atoms in the bottom-most 
layer and the hydrogen atoms terminating their dangling bonds, reaching a tolerance in the 
force of F max <1.0 mH/bohr. For the comparison of the energy advantages between the 
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FIG. 3: Energy advantage of germanium and silicon emitting structures compared with nonemitting 
structures as a function of the number of the oxygen atoms per unit cell, n. 



TABLE I: Energy advantage of germanium and silicon emitting structures compared with none- 
mitting structures in interface models, n indicates the number of the inserted oxygen atoms at the 
interface. All the values are in eV/unit cell. 





n = 3 


n = 6 


Quartz/Ge(100) 


-1.01 


-0.31 


Quartz/Si(100) 


-0.13 


3.94 



cases of germanium and silicon, the same models for Si/Si0 2 were also calculated. 

First, we examine the initial oxidation process of the substrates using the surface models. 
According to the oxidation process proposed by Kageshima and Shiraishi |7j, the silicon 
atom in the second layer will be emitted to release the stress. Figures [H(a) and []](b) show 
the examples for the case of n — 3, where n is the number of the inserted oxygen atoms. 
The hatched silicon atom is kicked out and goes elsewhere. When n is smaller than five, 
the dangling bonds remain at the interface after the emission. On the contrary, when six 
oxygen atoms are introduced, these dangling bonds are terminated by forming a Ge-O-Ge 
bond [see Fig. (He)]. Figure [3] shows the energy advantages of germanium and silicon as a 
function of the number of the inserted oxygen atoms. The energy advantage is denned as 
E n0 n{n) — (E emi {n) + where E non (n) and E emi (n) are the total energies of the nonemitting 
and emitted structure with n oxygen atoms inserted, respectively. In addition, \i is the 
chemical potential of the germanium or silicon in the bulk phase 18) . The energy advantage 



5 



TABLE II: Dispersion of bonding angles and thickness of oxidized region. The dispersions are 
calculated as X^ILi (109.5 — 6i) 2 /n, where 9i indicates the bonding angles around germanium (sil- 
icon) atoms, and n = 24 corresponding to the bonding angles of atoms in the first and second 
germanium (silicon) atomic layers per unit cell. The oxide region is the distance between the first 
and second germanium (silicon) atomic layers, ao indicates the experimental lattice constants of 
the bulk phase. 





Dispersion [degree 2 ] 


Thickness 


Ge 


283.4 


0.43a 


Si 


246.1 


0.40a 



of germanium for the atom emission is larger than that of silicon at n = 1 because the 
formation energy of monovacancy is significantly smaller in germanium (~ 1.9 eV) [3] 
than silicon (~ 3.3 eV) 20] . On the other hand, the energy advantage of silicon increases 
more rapidly as the oxygen atoms are inserted and become larger than that of germanium. 
The silicon-atom emission from the surface model is preferential when n is larger than 
three, while the germanium-atom emission is only at n — 6. We also examine the energy 
advantage of the atom emission from the interfaces as shown in Fig. [21 Table [I] summarizes 
the energy advantages for the emission from the interface, which are computed by the same 
manner with the case of the surface models. The energy advantage of germanium remains 
smaller than that of silicon. Although the quantitative discussion for the energy advantage 
is difficult because the chemical potentials of the atoms are simply computed from those in 
the bulk, our results indicate that the germanium atom preferentially stays in the substrate 
during the oxidation. In addition, if the germanium atoms are kicked out at n = 6, the 
two remaining dangling bonds are saturated by the reconstruction of the Ge-O-Ge bonding 
network as mentioned above. Since the event creating the dangling bonds at the interface 
scarcely happens, the defect density of the Ge/Ge02 interface becomes lower than that of 
the Si/Si0 2 interface. 

When oxygen atoms are sequentially inserted between Ge-Ge bonding networks from 
the surface, the formed oxide region has Ge-O-Ge bonding networks similar to that of the 
cristobalite of Ge02 bulk, and the bonding networks are attained lattice mismatch of 21 
%, which is almost the same as the lattice mismatch between Si and Si02- Therefore, the 
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difference in the lattice mismatch is not strongly related to the small energy advantage of 
the Ge/Ge02 interface. Then, we compare the dispersions of the bonding angles in the 
GeO"2 and SiC>2 regions at n = 4 which correspond to the 1-atomic-layer-oxidized surface. 
Table ITU shows the dispersions of the bonding angles from the tetrahedral structure of 109.5 
degree and the thickness of the oxidized region. The larger dispersion of the bonding angles 
around germanium means that the bonding angles around germanium are more flexible than 
that around silicon. Therefore, the Young modulus of GeO"2 is lower than that of Si02 in 
the good agreement with the experiment 2l|. Since the O-Si-0 bonding angles keep the 
tetrahedral structure rigidly, the thickness of the Si02 region barely increases to release the 
stress in the lattice mismatch while the Ge02 region elongates to the vacuum region. Thus, 
the high flexibility of the O-Ge-0 bonding networks causes the small energy advantages of 
the germanium-atom emission, and then leads to the lower defect density at the Ge/Ge02 
interface. Furthermore, the emitted atoms might cause the source of self-interstitials. The 
creation of these interstitials at the interface or in the substrate also becomes a serious 
problem because the interstitials cause the degradation of the MOS devices such as the 
dielectric leakage current and scattering centers of inversion carriers. Since the formation 
energy of the silicon self-interstitials (~ 3.5 eV) is almost equal to that of germanium self- 



interstitials, the 



the interface 
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arge energy advantage of silicon also gives rise to the self-interstitials at 



23j. Thus, our results imply that the Ge/Ge0 2 MOS field effect transistor 



exhibits the greater performance than the Si/Si02 one owing to the lower interface trap 
density at its interface as well as high carrier mobility. 

In conclusion, we have investigated the germanium- and silicon-atom emission during the 
oxidation process using the first-principles calculation. It is revealed that the germanium- 
atom emission hardly occurs compared with the silicon-atom emission because the high 
flexibility of the O-Ge-0 bonding networks decreases the stress due to the lattice mismatch 
at the Ge/Ge02 interface. Our study supports the experimental results that the lower 
interface trap density will be realizable in the Ge/Ge02 interface. Since we studied only 
limited models of surfaces and interfaces, not all of the details of the actual process are 
included. Nevertheless, we believe that our study will help to design the germanium based 
device as the beyond scaling devices. 
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